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We propose a systematic procedure to study a generalized linear sigma model which can give a 
physical picture of possible mixing between qq and qqqq low lying spin zero states. In the limit of zero 
quark masses, we derive the model independent results for the properties of the Nambu Goldstone 
pseudoscalar particles. For getting information on the scalars it is necessary to make a specific choice 
of terms. We impose two plausible physical criteria - the modeling of the axial anomaly and the 
suppression of effective vertices representing too many fermion lines - for limiting the large number 
of terms which are allowed on general grounds. We calculate the tree-level spectrum based on the 
leading terms in our approach and find that it prominently exhibits a very low mass isosinglet scalar 
state. Finally we point out that the low energy result for scattering of pions continues to hold in 
the general version of the model. 



Evidence has been accumulating [lj for a very light 
mass scalar-isoscalar particle, /o(600) as well as a pos- 
sible similar scalar- isospinor particle, 1<"q(800). As has 
been widely discussed, these may be joined with the well 
established scalars ao(980) and /o(980) to make a pu- 
tative light scalar nonet. The upside down mass order- 
ing of such a nonet suggests a four quark rather than a 
two quark structure; both qq-qq Q and qq-qq Q forms 
have been proposed. Either alternative would be of great 
importance for a full understanding of QCD in its non- 
perturbative low energy regime. The relation to the usual 
qq scalar mesons is of clear relevance in a such a picture. 
It has been suggested [ij that a mixing between the two 
quark and four quark nonets may help to better under- 
stand certain anomalies of the two quark nonet spectrum. 
The resulting picture is complicated and one may won- 
der, for example, whether the ordinary pions (believed 
to be of qq type) are chiral partners of the lighter four 
quark scalars or heavier two quark scalars. Historically, 
it has been of great value to study such questions in the 
framework of simple linear sigma models. Such a gener- 
alized linear sigma model was proposed in [5| and studied 
further in i and 0- These papers have suggested the 
plausibility of a situation in which the lightest, approxi- 
mate Nambu-Goldstone boson pseudoscalars are followed 
in ascending mass by scalars with relatively large four 
quark content. However, the model of interest may have 
many more terms than previously considered; for exam- 
ple if the interaction terms are restricted to be renor- 
malizable, there are 0] 21 chiral invariant terms and 21 
additional terms with the chiral transformation property 
of the QCD mass terms. In the present note we attempt 
to understand the essential structure more clearly, to dif- 
ferentiate between model dependent and model indepen- 
dent results as well as to suggest physical ways to choose 
the most important terms. As an aid we first simplify the 
analysis by setting the light quark masses to zero. It is 
accepted that this is a reasonable qualitative approxima- 
tion since the largest parts of the masses of all particles 



made of light quarks, other than the lightest 0~ octet, are 
expected to arise from spontaneous breakdown of chiral 
symmetry. 

The fields of our "toy" model consist of a 3 x 3 matrix 
chiral nonet field M , which represents qq type states as 
well as a 3 x 3 matrix chiral nonet field M', which rep- 
resents four quark type states. They have the decompo- 
sitions into scalar and pseudoscalar pieces: M — S + i(j), 
M' = S' + iqb' and behave under SU(3) L X SU(3) R trans- 
formations as M -> U L MUl and M' ->• U L M'U R . How- 
ever, the U(1)a transformation which acts at the quark 
level as q aL — > e m q aL , q a R — > e~" ' q a R distinguishes the 
two fields [B| according to 
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e 2w M, 



M' 



'M'. 



(1) 



Note that our treatment is based only on the symme- 
try structure and hence applies when M' is any linear 
combination of qq-qq and qq-qq type fields. We will be 
interested in the situation where non-zero vacuum values 
of S and S' may exist: (S£) = aS b a , (S'J>) = [35 b a , corre- 
sponding to an assumed SU(3)v invariant vacuum. The 
Lagrangian density which defines our model is 



C = - \tt (^M^Mt) - iTr (fyM'fyM *) 
- V (M,M')-V S B, 



(2) 



where Vo(M, M') stands for a general function made from 
SU(3)l x SU(3)r (but not necessarily U(1)a) invariants 
formed out of M and M'. The quantity Vsb which repre- 
sents the effective chiral symmetry breaking light quark 
mass terms will be set to zero here. 

The U(1)a transformation, which plays a special role 
in this model, suggests another useful simplification. In 
QCD there is a special instanton induced term- the '"t 
Hooft determinant" [U- which breaks the U(1)a sym- 
metry and can be modeled as det(M) + det(Aft). It 
thus may be natural to require all the terms to satisfy 
U(1)a invariance except for a particular subset which 
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could model the U(1)a anomaly. If one demands that, as 
in QCD, the U(1)a variation of the effective Lagrangian 
be proportional to the gluonic axial anomaly then a sim- 
ilar effective term £„ = -c 3 [ln(det(M)) - ln(det(Mt))P, 
where C3 is a numerical parameter, seems appropriate [9J. 

A lot of information concerning especially the pseu- 
doscalar particles in the model may be obtained in gen- 
eral without even specifying the terms in the potential. 
This may be achieved by studying "generating" equa- 
tions which arise from the demand that the infinitesimal 
symmetry transformations in the model hold. For the 
nine axial transformations one finds Q: 

to ||]+ - [S, ^]+ + (^,5)^(^0 = 
iroT fA' dVa o' dv o^ o n ( detM m (<i\ 



v detMt' 



To get constraints on the particle masses we will differ- 
entiate these equations once with respect to each of the 
two matrix fields <j> and <f>' and evaluate the equations 
in the ground state, taking into account the "minimum" 
conditions, (^) = and = 0. Further differen- 

tiations with respect to all four matrix fields will sim- 
ilarly yield "model independent" information on 3 and 
4 point vertices. We also require the Noether currents, 
(J° xial ) b a = ad„cf) b a + Pdpfil + ■ ■ ■, where the dots stand 
for terms bilinear in the fields. Using Eq.Q the squared 
mass matrix which mixes the degenerate two quark and 
degenerate four quark pseudoscalar octets is: 



{Ml) = to 



/3 2 /a 2 
-Pfa 



-p/a 
1 



(4) 



where = 



jrva\ 



Clearly, det(M^) = and the 



zero mass pion octet is a mixture of two quark and four 
quark fields. The transformation between the diagonal 
fields 7r + and tt' + and the original pion fields is defined 
as: 



Note that the physical higher mass pion state decouples 
from the axial current. Altogether the (eight) zero mass 
pseudoscalars are characterized by the three parameters 
a, (3 and y n . On the other hand, there are only two exper- 
imental inputs: — 131 MeV and the mass of 7r(1300), 
the presumed higher mass pion candidate. Thus, the in- 
teresting question of the four quark content of the pion 
has an inevitably model dependent answer. 

Next let us consider the "model independent" in- 
formation available for the two pseudoscalar SU(3) 
singlet states. This sector is related to the QCD 
axial anomaly. In the single M model, the anomaly 
can be modeled by the term C v mentioned above. 
In the M-M' model under consideration this form 
is no longer unique and it is natural to consider 
a generalization [10] in which ln(det(M)/det(M'l')) 
is replaced by 7i[ln(det(M)/det(M t ))] + (1 - 
7i)[ln(Tr(A/Af' t )/Tr(M'M t ))], where 71 is a di- 
mensionless parameter. Then the squared mass matrix 
which mixes the two SU(3) pseudoscalar singlet states is 
obtained as: 



(M 2 ) = 



8c 3 (27i+l)' 



3a- 



-z yo 



8c 3 (l- 



7 i)(2 7 i+1) 
3a/3 



-zoVo 4 

yo 



8c 3 (l- 



-71X271 + I) 
3a/3 
8c 3 (l-7l) 2 



-2/3/a and yo = (; 



d 2 V 



(8) 

Note that when 



Here z - ^ y/ ^ yu - \ g^g^, g , 
C3 is set to zero, making the entire Lagrangian U(1)a 
invariant, det(MQ) = 0. Then one of the singlet pseu- 
doscalars becomes, as well known, a Nambu Goldstone 
particle. This occurs in the large number of colors limit 
but we will not make that approximation here. 

In order to get information about the scalar meson 
masses and mixings as well as to complete the description 
of the pseudoscalars it is necessary to make a specific 
choice of interaction terms. To proceed in a systematic 
way we define the following quantity for each term, 



N = 2n + An', 



(9) 
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cos f 
sin f 



—sin 6 m 
cos Qtt 



which also defines the transformation matrix, 
explicit diagonalization yields: 



Rtt- 



tan( 







(5) 
The 



(6) 



which may be interpreted as the ratio of the four quark 
condensate to the two quark condensate in the underlying 
QCD. We see that the mixing between the two quark pion 
and the four quark pion would vanish if the four quark 
condensate were to vanish in this model. Rewriting the 
Noether current as {J^ lal )\ = F 7r d f ,Tr+ + F^d^i 
shows that 



F n , = 0. 



(7) 



where n and n' are respectively the number of M fields 
and the number of M' fields contained in that term. We 
shall restrict our choice to the lowest non-trivial value of 
N, which corresponds physically to the total number of 
quark and antiquark lines at each vertex. In addition to 
the two special terms which saturate the U(1)a anomaly 
already mentioned, this gives the leading (N—8) poten- 
tial 



V = - c 2 Tr(MM r ) + 
+ d 2 Tr (M'M't) 



clTr(MM^MM^) 

f e a 3 (e abc e def MSM b e M' f c 



h.c.) 

(10) 



where the dots stand for the U(1)a non-invariant 
terms. For simplicity, we have neglected the N=8 term, 
c|[Tr(M N'F)] 2 which is suppressed, in the single M 
model by the quark line rule. It may be noted that 
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the quantities det(M) and Tr(MM't) which enter into 
those two terms which saturate the U(1)a anomaly have 
iV=6. In this counting scheme, U(1)a invariant terms 
with N=12 (and higher) might be successively added to 
improve the approximation. The minimum equations for 
this potential are: 



dV 



2a[- 



c 2 



+ 4e«/3]=0, (11) 



0. 



(12) 



The U(1)a violating C3 terms do not contribute to these 
equations. Notice that a is an overall factor in Eq. (fTTj) 
so that, in addition to the physical spontaneous break- 
down solution where a^O there is a solution with a = 0. 
On the other hand, j3 is not an overall factor of Eq. ([12")) 
and it is easy to see that (3, which measures the "4 quark 
condensate" , is necessarily non-zero in the physical sit- 
uation where a is non-zero. From the specific form in 
Eq. (|10[) we find the mixing squared mass matrices for 
the degenerate octet scalars as well as the SU(3) singlet 
scalars: 



2 [-c 2 + 6cla 2 -2ef/3] 



2[" 



-4ae3 



c 2 + 6cf a' 
8ae% 



+ 4e§/3] 



-4ae§ 
2d, 



8ae§ 



2d, 



(13) 



(14) 



Now let us consider the comparison of the model with 
experiment. To start with there are 8 parameters (a, /3, 
C2, c?2, c|, eg, C3 and 71). The last two parameters ap- 
pear only in the mass matrices of the pseudoscalar SU(3) 
singlets and are conveniently discussed separately. The 
other six are effectively reduced to four by using the two 
minimum equations pip and (|12[) . As the correspond- 
ing four experimental inputs [1[ we take the non-strange 
quantities: 



m(0 + octet) 
m(0 + octct') 
m(0~octet') 



m[a (980)] = 984.7 ± 1.2 MeV 
m[a (1450)] = 1474 ± 19 MeV 
m[7r(1300)] = 1300 ± 100 MeV 
131 MeV (15) 



Evidently, the largest experimental uncertainty appears 
in the mass of 7r(1300); we shall consider the other masses 
as fixed at their central values and vary this mass in the 
indicated range. From studying the scalar SU(3) singlet 
states we find the consistency condition for positivity of 
the eigenvalues of their squared mass matrix, Eq. (TTJ 



The model predicts, as m[7r(1300)] varies from 1200 to 
1300 MeV, 



m(0 + singlet) 
m(0 + singlet') 



510 -> 28 MeV, 

1506 -> 1555 MeV. (17) 



to[tt(1300)] < 1302 MeV. 



(16) 



Clearly, the most dramatic feature is the very low mass 
of the lighter SU(3) singlet scalar meson. Of course, one 
expects the addition of quark mass type terms to modify 
the details somewhat. 

To calculate the masses of the SU(3) singlet pseu- 
doscalars we must diagonalize Eq.© with the specific 
choices of parameters yo = 2g?2 and zq = ^e^afdi corre- 
sponding to the potential of Eq. (TiTj)) . This enables us to 
fit in principle, for any choice of m[7r(1300)], the two pa- 
rameters C3 and 71 in terms of the experimental masses 
of T7(958) and one of the candidates ry(1295), 77(1405), 
?7(1475) and ?7(1760). However, it turns out that the posi- 
tivity of the eigenvalues of the matrix (Mq ) imposes addi- 
tional constraints on the choice of m[7r(1300)] in Eq. tfTB"]) . 
Furthermore the first two candidates for the heavier 77 are 
also ruled out on grounds of this positivity. For 77(1475) 
the allowed range of m[7r(1300)] is restricted to 1200 to 
1230 MeV. On the other hand, there is no additional re- 
striction if 77(1760) is chosen. If the choice of 77(1475) is 
made, the predicted range of m(0 + singlet) is narrowed 
from that given in Eq.lfT?)) to 510 -> 410 MeV. 

It is very interesting to see what the model has to 
say about the four quark percentages of the particles 
it describes. These percentages are displayed in Figfl] 
as functions of the precise value of the input parame- 
ter to[7t(1300)]. The pion four quark content (equal to 
100 sin 2 ^) is seen to be about 17 percent. Of course 
the heavier pion would have about an 83 percent four 
quark content. On the other hand, the octet scalar states 
present a reversed picture: the a (980) has a large four 
quark content while the ao(1450) has a smaller four quark 
content. The very light and the rather heavy + singlets 
are about maximally mixed, having roughly equal con- 
tributions from the 4 quark and 2 quark components. 

The perhaps more plausible scenario in the case of the 
0~ singlets takes ?7(1475) as the heavy 0~ singlet state. 
Fig. Q] shows that there are two solutions for each value 
of m[7r(1300)]; the dotted line gives a mainly qq content 
while the solid line gives a mainly four quark content. 
Note that this scenario does not allow to[7t(1300)] to be 
higher than about 1230 MeV. The choice (not shown) of 
77(1760) as the partner of 77(958) also leads to two solu- 
tions with small and large two quark content. 

There are two reasons for next briefly discussing the 
pi-pi scattering in this model. First, since the iso-singlet 
scalar resonances above are being considered at tree level, 
one expects, as can be seen in the single M model also dis- 
cussed in [5j and at the two flavor level in 11 1 , that uni- 
tarity corrections for the scattering amplitudes will alter 
their masses and widths. Second, since the pion looks un- 
conventional in this model (having a non-negligible four 
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1.24 1.26 
m [n(1300)] (GeV) 



that while the (lighter) pion is mainly two quark, the 
lighter scalars have very large four quark components. 
This is perhaps the opposite of what one might initially 
think and is related to the characteristic mixing pattern 
emerging in a transparent form here. 
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FIG. 1: Four quark percentages of the pion (dashed line), 
the do (980) (top long-dashed line), the very light + singlet 
(dotted-dashed line) and the 77(958) in the scenario where the 
higher state is identified as the 77(1475) (curve containing both 
solid and dotted pieces) as functions of the undetermined in- 
put parameter, m[7r(1300)]. Note that there are two solutions 
for the 77(958): the dotted curve choice gives it a predominant 
two quark structure and the solid curve choice, a larger four 
quark content. 



quark component) one might worry that the fairly pre- 
cise "current algebra" formula for the near to threshold 
scattering amplitude might acquire unacceptably large 
corrections. In the present massless pion model, this for- 
mula [12| should read, 



A(s,t,u) = 2s/Fl 



(18) 



where A(s, t, 11) is the conventional amplitude term ex- 
pressed in terms of the Mandelstam variables s, t and u. 
To obtain A(s, t, u) in the present model, one needs the 
four point vertices involving the pseudoscalar octet fields 
as well as the three point vertices involving two pseu- 
doscalar octet fields and one scalar field. It turns out 
fl3l ] that the result Eq. ([l~8"]) follows in a "model indepen- 
dent" way just by using the generating Eq.Q: the four 
point vertices can be related to the three point vertices, 
which can in turn be related to the two point vertices 
(masses). For example, the three point vertices involving 
the SU(3) singlet scalars can be related to the scalar and 
pseudoscalar squared mass matrices as: 



/ d 3 V 

^d(4>'i) A d(4' 1 2 )Bd(S )H 



AH 



AH- 



(19) 



Here the capital Latin subscripts refer to summation over 
the unprimed and primed fields, (M%) is given in Eq.Jl]) 
and (Xq) is the model independent version of Eq. (fT4)l . 
It is interesting to note that the current algebra theo- 
rem will "tolerate" any amount of four quark component 
in the massless pion. The present model clearly shows 
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